Current polygenic scores (PGS) have limited predictive power for fracture risk. To improve genetic prediction, we developed and evaluated a novel metaPGS combining genetic information from multiple fracture-related traits. We derived individual PGS from genome-wide association studies of 16 fracture-related traits and employed an elastic-net logistic regression model to examine the association between the 16 PGSs and fractures. An optimal metaPGS was constructed by combining 11 significant individual PGSs selected by the elastic regularized regression model. We evaluated the predictive power of the metaPGS alone and in combination with clinical risk factors recommended by guidelines. The discrimination ability of metaPGS was assessed using the concordance index. Reclassification was assessed using net reclassification improvement (NRI) and integrated discrimination improvement (IDI). The metaPGS had a significant association with incident fractures (HR 1.21, 95% CI 1.18-1.25 per standard deviation of metaPGS), which was stronger than previously developed bone mineral density (BMD)-related individual PGSs. Models with PGS_FNBMD, PGS_TBBMD, and metaPGS had slightly higher but statistically non-significant c-index than the base model (0.640, 0.644, 0.644 vs. 0.638). However, the reclassification analysis showed that compared to the base model, the model with metaPGS improves the reclassification of fracture. The metaPGS is a promising approach for stratifying fracture risk in the European population, improving fracture risk prediction by combining genetic information from multiple fracture-related traits.